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Abstract 

We consider a finite quantum system S coupled to two environments of different 
nature. One is a heat reservoir 1Z (continuous interaction) and the other one is a 
chain C of independent quantum systems £ (repeated interaction). The interactions 
of S with 1Z and C lead to two simultaneous dynamical processes. We show that for 
generic such systems, any initial state approaches an asymptotic state in the limit 
of large times. We express the latter in terms of the resonance data of a reduced 
propagator of S + 1Z and show that it satisfies a second law of thermodynamics. 
We analyze a model where both S and £ are two-level systems and obtain the 
asymptotic state explicitly (lowest order in the interaction strength) . Even though 
7Z and C are not direcly coupled, wc show that they exchange energy, and we find 
the dependence of this exchange in terms of the thermodynamic parameters. 

We formulate the problem in the framework of W* -dynamical systems and 
base the analysis on a combination of spectral deformation methods and repeated 
interaction model techniques. Wc do not use master equation approximations. 

1 Introduction 

Over the last years, the rigorous study of equilibrium and non-equilibrium quantum 
systems has received much and renewed attention. While this topic of fundamental 
interest has a long tradition in physics and mathematics, conventionally explored via 
master equations [9, 6], dynamical semi-groups [3, 6] and algebraic scattering theory 
[32, 15], many recent works focus on a quantum resonance theory approach. The latter 
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has been applied successfully to systems close to equilibrium [17, 26, 27, 28] and far 
from equilibrium [18, 25]. In both situations, one of the main questions is the (time-) 
asymptotic behaviour of a quantum system consisting of a subsystem S interacting 
with one or several other subsystems, given by thermal reservoirs Hi, . . . ,TZ n . It has 
been shown that if S + 1Z starts in a state in which the reservoir is in a thermal state 
at temperature T > far away from the system S, then S + 1Z converges to the joint 
equilibrium state at temperature T, as time t — > oo. This phenomenon is called re- 
turn to equilibrium. (See also [22] for the situation where several equilibrium states 
at a fixed temperature coexist.) In case S is in contact with several reservoirs having 
different temperatures (or different other macroscopic properties), the whole system 
converges to a non- equilibrium stationary state (NESS). The success of the resonance 
approach is measured not only by the fact that the above-mentioned phenomena can 
be described rigorously and quantitatively (convergence rates), but also by that the 
asymptotic states can be constructed (via perturbation theory in the interaction) and 
their physical and mathematical structure can be examined explicitly (entropy pro- 
duction, heat- and matter fluxes). One of the main advantages of this method over 
the usual master equation approach (and the related van Hove limit) is that it gives a 
perturbation theory of the dynamics which is uniform in time t > 0. While the initial 
motivation for the development of the dynamical resonance theory was the investiga- 
tion of the time-asymptotics, the method is becoming increasingly refined. It has been 
extended to give a precise picture of the dynamics of open quantum systems for all 
times t > 0, with applications to the phenomena of decoherence, disentanglement, and 
their relation to thermalization [27, 26, 28, 23]. An extension to systems with rather 
arbitrary time-dependent Hamiltonians has been presented in [29] (see also [2] for time- 
periodic systems). A further direction of development is a quantum theory of linear 
response and of fluctuations [20]. 

In certain physical setups, the reservoir has a structure of a chain of independent 
elements, C = 8\ + £ % + ■ • • . An example of such a system is the so-called "one-atom 
maser" [24] , where S describes the modes of the electromagnetic field in a cavity, inter- 
acting with a beam C of atoms £j, shot one by one into the cavity and interacting for a 
duration tj > with it. A mathematical treatment of the one-atom maser is provided 
in [14]. Another instance of the use of such systems is the construction of reservoirs 
made of "quantum noises" by means of adequate scaling limits of the characteristics of 
the chain C and its coupling with S, which lead to certain types of master equations 
[1, 7, 6, 4, 5]. The central feature of such systems is that S interacts successively with 
independent elements £j constituting a reservoir. This independence implies a marko- 
vian property which simplifies the mathematical treatment considerably. In essence it 
enables one to express the dynamics of S at time t = t\ + • • • + tn by a propagator of 
product form M\{t\) ■ ■ ■ Mn{tn)-> where each Mj(rj) encodes the dynamics of S with 
a fixed element £j. In case each element £j is physically the same and each interaction 
is governed by a fixed duration r (and a fixed interaction operator), the dynamics is 
given by M(t) n and the asymptotics is encoded in the spectrum of the reduced dy- 
namics operator M(r) [11]. An analysis for non-constant interactions is more involved. 
It has been carried out in [12, 13] for systems with random characteristics (e.g. ran- 
dom interaction times). See also [30] for related issues. In both the deterministic and 
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the random settings, the system approaches a limit state as t — > oo, called a repeated 
interaction asymptotic state (RIAS), whose physical and mathematical properties have 
been investigated explicitly. 

In the present work we make the synthesis of the above two situations. We consider 
a system S interacting with two environments of distinct nature. On the one hand, 
S is coupled in the repeated interaction way to a chain C = £ + £ + ••• , and on the 
other hand, S is in continuous contact with a heat reservoir TZ. It is assumed that C 
and TZ do not interact directly. Our goal is to construct the asymptotic state of the 
system and to analyze its physical properties. In particular, we present in Section 1.2 
our results on the convergence to, and form of the asymptotic state, in Section 1.3 
the thermodynamic properties of it, and in Section 1.4 we present the analysis of an 
explicit model. 

1.1 Description of the system 

The following is a unified description of S, TZ, C in the language of algebraic quantum 
statistical mechanics (we refer the reader to e.g. [31] for a more detailed exposition). 
The Hilbert spaces of states of each of the subsystems # = S,1Z,£ are Tijt. The 
respective observables form von Neumann algebras 9Jt# C B(Tt#). We assume that 
dimHs < oo and dimTC^ may be finite or infinite. TZ being a reservoir, its Hilbert 
space is assumed to be infinite-dimensional, dim?^ = oo. The free dynamics of each 
constituent is generated by Liouville operators L#, i.e., the Heisenberg evolution of 
an observable A £ 9Jt# at time t is given by e itL # Ae~ ltL # . In each Hilbert space we 
pick a normalized reference state which determines the macroscopic properties of 
the systems. 1 These reference vectors are invariant, e~ ltL # i &# = and they are 
cyclic and separating for 9Jt# [10]. Typically, the ty# are chosen to be the equilibrium 
states at any fixed temperature T# > 0. The Hilbert space Tic of the chain is the 
infinite tensor product of factors Tie, taken with respect to the stabilizing sequence 

In summary, the non-interacting system is given by a Hilbert space 

Ti = H s <8> H n <g> H c 
and its dynamics is generated by the Liouvillian 

L = L s + L n + Y,L £k . (1.1) 

k>i 

Here we understand that Lg k acts as the fixed operator Lg on the /c-th factor of Tic, 
and we do not display obvious factors 11. 

The operators governing the couplings between S and £ and S and TZ are given by 

V S£ £Wl s ®W £ and V&r € Wis ® SDfo 

1 In other words, it determines the folium of normal states. If the Hilbert space is finite-dimensional 
then the set of normal states is unique, but for infinite systems different classes of normal states are 
determined by different macroscopic parameters, such as the temperature. 
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respectively, and the total interaction is 

V(\) = \ 1 Vsii + *2V S eeWls®<m n ®Wl s , (1.2) 
where Ai, A2 are coupling constants (A = (Ai, A2)). The full (Schrodinger) dynamics is 

ip^U(m)ip, (1.3) 

where U(m) is the unitary map 

JJ( m \ — e -iT(io+V r m) e -ir(-L +V r m _i) . . . g-ir^o+^i)^ ^.4) 

r > being the time-scale of the repeated interaction and Vk being the operator V(A), 
(1.2), acting nontrivially on Tig, Tin and the fc-th factor Tig of Tic ( we will also write 
L m = Lq + V m ). We discuss here the dynamics (1.3) at discrete time steps mr only, 
a discussion for arbitrary continuous times follows in a straightforward manner by 
decomposing t = mr + s, s 6 [0,r), see [11]. 

Explicit form of finite systems and thermal reservoirs. 

(A) Finite systems. We take 5 (and possibly £ ) to be finite. The Hamiltonian 
of 5 is given by H$, acting on ()£. The (Gelfand-Naimark-Segal) Hilbert space, the 
observable algebra and the Liouville operator are given by 

Hs = t)s®*)s, fms = B(n s )<si, L s = H S ®1-1® H s . 
The reference state is chosen to be the trace state, represented by 

^ dim'Hs 

\/dimHs ~[ 

where {<Pj} is an orthonormal basis of i)s diagonalizing H$- 

(B) Thermal reservoirs. We take 1Z (and possibly £) to be a thermal reservoir of free 
Fermi particles at a temperature T > 0, in the thermodynamic limit. Its description 
was originally given in the work by Araki and Wyss [8]; see also [17] and [29], Appendix 
A, for an exposition close to ours. The Hilbert space is the anti-symmetric Fock space 

Wtt = r_«,) :=0P-[L 2 ft)]^ 

n>0 

over the one-particle space 

f) = L 2 (]R,©), (1.5) 

where V- is the orthonormal projection onto the subspace of anti-symmetric functions, 
and <S is an 'auxiliary space' (typically an angular part like L 2 (S 2 )). In this represen- 
tation, the one-particle Hamiltonian h is the operator of multiplication by the radial 
variable (extended to negative values ) s G M of (1.5), 

h = s. 
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The Liouville operator is the second quantization of h, 

n 

L n = dT(h):=@Y, h 3> (1-6) 

n>0 j=l 

where hj is understood to act as h on the j-th factor of V- [L 2 ^)]®^ 1 and trivially on 
the other ones. 

The von Neumann algebra Win is the subalgebra o£B(?in) generated by the thermal 
fermionic field operators (at inverse temperature /?), represented on Tin by 

<p(0p) = -^[a*(9p) + a,{gp)]. 

Here, we define for g G L 2 (R + , (5) 



9p(s) 



g(s) if s > 

e -f3s + 1 1 g(- S ) if s < 0. 



We choose the reference state to be thermal equilibrium state, represented by the 
vacuum vector of Tin, 

We provide a precise derivation of the above formalism starting from the usual descrip- 
tion of a reservoir of non-interacting and non-relativistic fermions in Section 1.4. 



1.2 Convergence to asymptotic state 

One of our main interests is the behaviour of averages p{U (m)* O m U (m)) as m — ► oo, 
where p is any (normal) initial state of the total system, and where O m is a so-called 
instantaneous observable [11, 12, 13]. Let Asu £ 9Jl$ ® 9Jtn and let Bj G Wig, j = 
—I, . . . , r, where l,r > are integers. The associated instantaneous observable is 

O m = A sn ^ (Bj_ m ) eTl (1.7) 

where "&j{B) is the observable of SOT which acts as B on the j-th. factor of Tic, and 
trivially everywhere else ($j is the translation to the j-th factor). An instantaneous 
observable is a time-dependent one. It may be viewed as a train of fixed observables 
moving with time along the chain C so that at time m it is "centered" at the m-th factor 
Tis of Tic j on which it acts as Bq . If O acts trivially on the elements of the chain, then 
the corresponding instantaneous observable is constant and O m = O. However, in 
order to be able to reveal interesting physical properties of the system, instantaneous 
observables are needed. For instance observables measuring fluxes of physical quantities 
(like energy, entropy) between S and the chain involve instantaneous observables acting 
non-trivially on 8 m and on £ m +ii which corresponds to nontrivial A$-ji and Bo, B\. 
We denote the Heisenberg dynamics of observables by (see (1.4)) 

a m {O m ) = U(m)*O m U(m). (1.8) 
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The total reference vector 

*0 = *5 ® ^TZ ® (1-9) 

is cyclic and separating for the von Neumann algebra 9JT. We also introduce, for later 
purposes, the projector Psn = ^sn ® l^cK^cli which range we often identify with 
'H-sil = T~(-s %'Hti- Let J and A be the modular conjugation and the modular operator 
associated to the pair (971, ^o) [10]. In order to represent the dynamics in a convenient 
way (using a so-called C-Liouville operator), we make the following assumption. 

HI The interaction operator V(X), (1.2), satisfies A 1/2 y(A)A" 1/2 G M s (g> 9% <g> Wl £ . 

Since we will be using analytic spectral deformation methods on the factor TC-ji of 
7i, we need to make a regularity assumption on the interaction. Let R 3 9 i— ► T(#) S 
B(7in) be the unitary group defined by 

T(#) =r(e- eas ) on T_(L 2 (M,0)), (1.10) 

where for any / G L 2 (R, U5), 

(e-«>V)(s) = /(*-*). 

In the following, we will use the notation 

T(0) = t s ® T{9) (8) 1 £ 

for simplicity. Note that T{9) commutes with all observables acting trivially on TC-jz, in 
particular with Psn- Also, we have T{9)^n = ^>ti for all 9. The spectral deformation 
technique relies on making the parameter 9 complex. 

H2 The coupling operator Wsn '■= Vsn — JA^Vs-^A -1 / 2 J is translation analytic 
in a strip Kg = {z : < < 9q} and strongly continuous on the real axis. 
More precisely, there is a 9q > such that the map 

R B e^ T- 1 (9)w sn T(9) = w sn (9) em s ® m n , 

admits an analytic continuation into 9 £ Kg which is strongly continuous as 
Q9 I 0, and which satisfies 

sup \\Wsn(e)\\ < oo. 
o<se<e 

Let O m be an instantaneous observable (1.7). We say that O m is an analytic observable 
if 

r(fl)- 1 o m * (i-ii) 

has an analytic extension to 9 £ ng which is continuous on the real axis. Evidently, 
since T acts on TL-ji only, this is equivalent with T{9)~ 1 Asti^o having such an extension. 

Finally we present a 'Fermi golden rule condition' which guarantees that the sub- 
systems are well coupled so that the physical phenomena studied are visible at lowest 
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nontrivial order in the pertrubation Ai, A2. This is a very common hypothesis which is 
most often verified in concrete applications. To state it, we mention that the evolution 
is generated by so-called reduced dynamics operators [11] M(A) acting on the reduced 
space TLsii (where the degrees of freedom of C have been 'traced out'). In this paper, 
we analyze the spectrally deformed operators 

M e {\) =T{dy l M{\)T(6). 

We show in Corollary 2.3 that Mg{X) has an analytic extension into the strip Kg , in 
the sense of H2 above. We show that 1 is an eigenvalue of M${X) for all 6, A and 
that, for small couplings A = (Ai, A2) , the spectrum of Mg(X) must lie in the closed 
unit disk. The latter fact is true because Mg(X) is the analytically translated (one- 
step) propagator of a reduced unitary dynamics, although we can only prove it in a 
perturbative regime. The former fact can be seen as a normalization (the trace of the 
reduced density matrix of S + 1Z equals unity at all times). Since the propagator at 
time step m is represented by a power of Mg(X), it is not surprising that convergence 
to a final state is related to the peripheral eigenvalues of Mg(X). The following Fermi 
golden rule condition is an ergodicity condition ensuring the existence of a unique limit 
state. 

FGR There is a 9\ E k$ and a Ao > (depending on Q\ in general) such that for all 
A with < |A| < Ao, a(Mg 1 (X)) (spectrum) lies inside the complex unit disk, 
and a(Mg 1 (X)) n S = {1}, the eigenvalue 1 being simple and isolated. (§ is the 
complex unit circle.) 

This condition is verified in practice by perturbation theory (small A). It is also 
possible to prove that if the spectral radius of Mg 1 (X) is determined by discrete eigen- 
values only, then the spectrum of Mg 1 (X) is automatically inside the unit disk (see 
Proposition A. 3). Since the spectrum is a closed set the FGR condition implies that 
apart from the eigenvalue 1 the spectrum is contained in a disk of radius e~ 7 < 1. 

Theorem 1.1 (Convergence to asymptotic state) Assume that assumptions HI, 
H2 and FGR are satisfied. Then there is a Ao > s.t. if < \X\ < Ao, the following 
holds. Let p be any normal initial state on 9Jt, and let O m be an analytic instantaneous 
observable of the form (1.7). Then 

r 

m lim o ^(a m (O m )) = p+,x(Psr a l+1 {A sn ®° = _, B^Psr) \\^e\B^ £ ), (1.12) 

m ™ > °° j= i 

where p+ t \ is a state on Wis ® ^Tl> Psiz is the orthogonal projection onto the subspace 
Tls®'ri / R> an d where a 1 is the dynamics (1.8). Moreover, for analytic A G 9Jt ; we have 
the representation 

P+, x (PskAPsk) = (r ei ( x )\ T (^y lp snAP S Tz^s ® (1.13) 

where ipg (X) is the unique invariant vector of the adjoint operator [Mg 1 (X)]* , normal- 
ized as (^ i (A)|^'o) = I- 
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Remark. The operators Bj with j > 1 measure quantities on elements £ m +j which, at 
time m, have not yet interacted with the system S. Therefore they evolve independently 
simply under the evolution of S m+ j. For large times m — > oo, the elements of the chain 
approach the reference state (because the initial state is normal), and the latter is 
stationary w.r.t. the uncoupled evolution. This explains the factorization in (1.12). 
As a special case of Theorem 1.1 we obtain the reduced evolution of S + 1Z. 

Corollary 1.2 Assume the setting of Theorem 1.1. Then 

lim p(a m (A S n)) =P+,x{Asn)- 

1.3 Thermodynamic properties of asymptotic state 

The total energy of the system is not defined, since 1Z and C are reservoirs (and typi- 
cally have infinite total energy). However, the energy variation is well defined. More 
precisely, the formal expression for the energy at time m, a m (Lo + $m(V)) (see (1.1), 
(1.2)), undergoes a jump 

AE to \m) = a m+1 (L + $ m+1 (V))-a m (L + $ m (V)) 

= \2a m ($m + l(V S £)-$m(V S £)), 

as time passes the moment mr. The variation AE tot (m) is thus an instantaneous 
observable. In applications this observable is analytic and hence we obtain under the 
conditions of Theorem 1.1 (see also [11]) that 

dE™:= lim -p(AE tot (m)) = p + ^jT), 

m— +oo Tfl 

where j+ l = V — ct T (V) is the total energy flux observable. The quantity dE^ 1 repre- 
sents the asymptotic energy change per unit time r of the entire system. In the same 
way we define the variation of energy within the system S, the reservoir 1Z and the 
chain C between times m and m + 1 by 

AE s (m) 
AE n (m) 
AE c {m) 

These variations can be expressed in terms of commutators [Vss,L#] and [Vsn, 
where # = S,£,TZ. Since [Vse, L#] acts on S + £ only, it is an analytic observable (see 
sentence after (1.11)). We make the following Assumption. 

H3 The commutators [Vsk,L#], where # = S,£,1Z, are analytic observables in Wl. 

We can thus apply Theorem 1.1 to obtain (see Section 4) 

dE* := lim -p(AE*(m)) = p +tX (j#), # = S,TZ,C, (1.14) 



= a m+1 (L s )-a m (L s ), 

= a m+1 (L n )-a m (L n ), 

= a m+ \L £m+1 )-a m {L £m+1 ). 
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where j& are explicit 'flux observables' (c.f. (4.3)-(4.5)). We show in Proposition 4.1 
that j^* = j+ + j+ + and that p+ = 0. It follows immediately that 

dE^ = dE% + dE%. (1.15) 

The total energy variation is thus the sum of the variations in the energy of C and 1Z. 
The details of how the energy variations are shared between the subsystems depends 
on the particulars of the model considered; see below for an explicit example. 



Next we consider the entropy production. Given two normal states p and po on 9JT, 
the relative entropy of p with respect to po is denoted by Ent(p|po)- (This definition 
coincides with the one in [11] and differs from certain other works by a sign; here 
Ent(p|p ) > 0). 

We examine the change of relative entropy of the state of the system as time evolves, 
relative to the reference state po represented by the reference vector ^o, see (1.9). For 
a thermodynamic interpretation of the entropy, we take the vectors # = S, £, 1Z 
to represent equilibrium states of respective temperatures (3$, fls, fin- We analyze the 
change of relative entropy 

AS(m) = Ent(p o a m \p ) - Ent(p|p ) 

proceeding as in [11]. We show in Section 4 (see (4.6)) that 

dS+ : = hm = {(3 n - (3 £ )dE^ + fedE™. 

Combining this result with (1.15), we arrive at 

Corollary 1.3 The system satisfies the following asymptotic 2nd law of thermody- 
namics, 

dS+ = Pedrf + foidE*. 



1.4 An explicit example 

We consider S and £ to be two-level systems. The observable algebra for S and for £ 
is %s = 2lg = M 2 (C). Let E s , Eg > be the "excited" energy level of S and of £, 
respectively. Accordingly, the Hamiltonians are given by 

HS = ( Es) k£ = ( Be)' 

The dynamics are ol s {A) = e lth sj{ e -^hs an( j Q,|(yl) = e lth£ Ae~ lthe . We choose the 
reference state of £ to be the Gibbs state at inverse temperature /?£, i.e. 

PPe,eW = - L y S where Z PejE = Tr( e -^), 
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and we choose (for computational convenience) the reference state for S to be the 
tracial state, ps(A) = ^Tr(A). The interaction operator between S and an element £ 
of the chain is defined by \2Vss-, where A2 is a coupling constant, and 

vse ■= as <8> a* £ + a* s <g) as ■ 

The above creation and annihilation operators are represented by the matrices 

a * = ( I ) and a # = ( \ ) • 

To get a Hilbert space description of the system, one performs the Gelfand-Naimark- 
Segal (GNS) construction of (%Ls,Ps) and (2lf, Pp e ,s)i see e -g- [10, 11]. In this repre- 
sentation, the Hilbert spaces are given by 

H s = Us = C 2 ® C 2 , 

the Von Neumann algebras by 

Ms = Ms = M 2 (C) ® i C 2 c B(C 2 ® C 2 ), 
and the vectors representing ps and pp £ ,s are 

* 5 = (|0) ® |0> + |1) ® |1» , *g = (|0) 8) |0) + e"^/ 2 |l) ® |1)) . 

In other words, = (^5, (-4 ® and pp £ ,s(A) = (tps, (A ® H)#). Above, |0) 

(resp. |1)) denotes the ground (resp. excited) state of hg and hs- For shortness, in the 
following we will denote for \i) <g> i, j = 0,1. The free Liouvilleans Ls and Ls 
are given by 

L s = h s ® 1 C 2 - 1c 2 ® h s , L £ = h e ® 1 C 2 - 1 C 2 ® /ig 
and the interaction operator Vss is 

Vse = {a s ® i c2 ) ® (a| ® l c2 ) + (aj ® l c2 ) ® (ag ® l c2 ). 

For the reservoir, we consider a bath of non-interacting and non-relativistic fermions. 
The one particle space is fyn = L 2 (M, 3 ,d 3 k) and the one-particle energy operator h-jz is 
the multiplication operator by \k\ 2 . The Hilbert space for the reservoir is thus r_(fj^) 
and the Hamiltonian is the second quantization dT(hn) of h-ji (see (1.6)). The algebra of 
observables is the C*-algebra of operators 21 generated by {a*(/) | / £ rj-^.} where a/a* 
denote the usual annihilation/creation operators on T_(i)Ti). The dynamics is given by 
Tf (a (/)) = a&(e lth f), where h is the Hamiltonian of a single particle, acting on f). It 
is well known (see e.g. [10]) that for any f}ji > there is a unique (rf, (3)— KMS state 
/9g K on 21, determined by the two point function pp n (a*(f)a(f)) = (f, (1 + e^ R - h ' R -)~ 1 f) , 
and which we choose to be the reference state of the reservoir. Finally, the interac- 
tion between the small system S and the reservoir is chosen of electric dipole type, 
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i.e. of the form vsn = ( a s + a-s) ® <Pft(f) where / E rj^ is a form factor and 

<p{J) = jiW) + «*(/))■ 

We know explain how to get a descrpition of the reservoir similar to the one given 
in Section 1.1. As for S and £, the first point is to perform the GNS representation 
of (21, P/3 n ), so called Araki-Wyss representation [8]. Namely, if Q denotes the Fock 
vacuum and N the number operator of r_(rj-fc), the Hilbert space is given by 

Ti n = r_(L 2 (M 3 ,d 3 A;)) ®r_(L 2 (R 3 ,d 3 £;)), 

the Von-Neumann algebra of observables is 

SWW = VT/3 (21)" 

where 

*„(«(/)) = a (^/) ® 1 + (-U* ® (^/) =: 
^(a*(/)) = a* (^==/) ® 1 + (-1)^ « a (^/) == «$(/), 
the reference vector is 

and the Liouvillean is 

= dT(h n ) dT{h n ). 

We then consider the isomorphism between L 2 (]R 3 ,d 3 /c) and L 2 (IR + x S 2 , -^-drdcr) ~ 

L 2 (IR + , -^dr; 0), where = L 2 (5 2 ,d<r), so that the operator hn (the multiplication 
by \k\ 2 ) becomes multiplication by r S ]R + (i.e. we have r = |fc| 2 ). The Hilbert space 
Ti-ji is thus isomorphic to 

L 2 (M + , ^dr; 0)J ® T„ ^L 2 (R+, ^dr; 0)J . (1.16) 

Next we make use of the maps 

a #(/) ®l^ a # (/©0), (-1)* (8) a # (/) ^ a # (0 /) 
to define an isometric isomorphism between (1.16) and 

L 2 (R + , ^-dr; 0) © L 2 (R+, ^dr; 0) 

A last isometric isomorphism between the above Hilbert space and 

Tin :=T_ (L 2 (M,ds;0)) 

is induced by the following isomorphism between the one-particle spaces L 2 (M+,^dr;0) 
L 2 (E+, ^dr; 0) and L 2 (R, ds; 0) =: fj 

|s| 1/4 / f(s) if s > 0, 



f ® g >—> h, where /i(s) 



^2 I g(s) ifs < 0. 
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Using the above isomorphisms, one gets a descrition of the form given in Section 
1.1 for the reservoir 1Z. In this representation, the interaction operator Dsn becomes 

Vsn = (<r x ® W ® vUp-r) e^s® Win, 

where a x = as + is the Pauli matrix and fp n G t) = L 2 (M, ds; L 2 (S 2 , da)) is related 
to the initial form factor / £ L 2 (R 3 ,d 3 fe) as follows 

(fT) = 71 Vl + e^ 1 f Ws a) if S < 0. ( L1? ) 

As mentioned at the beginning of the introduction the situation where 5 is inter- 
acting with 7Z or C alone has been treated in previous works [17, 29] and [11]. If S is 
coupled to 1Z alone, then a normal initial state approaches the joint equilibrium state, 
i.e. the equilibrium state of the coupled system S + 1Z at temperture (5^ , with speed 
e -mr7th ( we consider discrete moments in time, t = tot to compare with the repeated 
interaction situation). If S is coupled to C alone, initial normal states approach a re- 
peated interaction asymptotic state, which turns out to be the equilibrium state of 5 
at inverse temperature /3' £ where 

0e = /fe|r. ( L18 ) 
and with speed e _mT7ri . The convergence rates are given by 



7th = Kl^+0{\\), with 7 w = -y/Es\\f(VEs)\\v (1-19) 

7n = \bS ] +0(\*), with 7 W = rsinc 2 ( l&LzM ) , (1 . 20 ) 

where sinc(x) = sin(a;)/x and ||/(\/^s)lll : = / \f(\/Es o")| 2 do\ 

75 2 



'5 2 

In order to satisfy the translation analyticity requirement H2, we need to make 
some assumption on the form factor /. Let 1(5) = {z £ C, |3 : (^)| < 5}. We denote by 
H 2 (5) the Hardy class of analytic functions h : 1(5) — > <5 which satisfy 



H*(S) ••= SU P / \\h(s + i9)\\%ds < oo. 
W\<5 JR 



\9\<S 

H4 Let f be defined by (1.17), with fa = 0. There is a 5 > s.t. e-^ s / 2 / (s) € 
H 2 (5). 

Proposition 1.4 (Asymptotic state of S) Assume f satisfies H4, ||/(\AEs)||<g ^ 
anc/ t(Ez — .E^) ^ 27rZ*. Then the asymptotic state p +t \ is given by 

P+,x = {iPp^s + (1 - 7)pp' s ,s) ® Pfe,^ + 0(A), 

where pp # is the Gibbs state of#, # = S,1Z, at inverse temperature (3 and where 7 is 
given by 

\2 /v (2) , \2,J2) ' 
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Remark. The fact that the asymptotic state p+ s \ is a convex combination of the two 
asymptotic states P/3 n) s and P/3' g ,s holds only because the system S is a two-level system 
and is not true in general. 

Using (4.2)-(4.5), Corollary 1.3 and Proposition 1.4, an explicit calculation of the 
energy fluxes and the entropy production for this concrete model reveals the following 
result. 

Proposition 1.5 Assume that \\f(y/Es)\\ e / and t(E £ - E s ) £ 2irZ* . Then 
dE c + = KEe(e-P* E s-e-K E A+0(\ 3 ), 
dE* = Kj E lS (e-^-e-^)+0(A 3 ), 
dE'f = k(E £ - E s ) (e-^s _ e -P'e E s\ + o(A 3 ), 
dS + = k(P' b E s - PkE s ) (e~^ Es - e"^) + 0(A 3 ), 



where 



\2 (2) ,2 (2) 
7 -l 7 -l A l7th A 27 ri 



j Pk,S"P £! S 2 (2) a2 (2)" 
A l7th + A 27ri 

Remarks. 1. The constant k is positive and of order A 2 . Moreover it is zero if at least 
one of the two coupling constants vanishes (we are then in an equilibrium situation and 
there is no energy flux neither entropy production). 

2. The energy flux dE+ is positive (energy flows into chain) if and only if the 
reservoir temperature Tn = is greater than the renormalized temperature Tg = 
(i'f 1 of the chain, i.e. if and only if the reservoir is "hotter". A similar statement holds 
for the energy flux d-E^ of the reservoir. Note that it is not the temperature of the 
chain which plays a role but its renormalized value (1.18). 

3. When both the reservoir and the chain are coupled to the system S (A1A2 7^ 0) 
the entropy production vanishes (at the main order) if and only if the two temperatures 
T-ji and T£ are equal, i.e. if and only if we are in an equilibrium situation. Once again, 
it is not the initial temperature of the chain which plays a role but the renormalized 
one. 

4. The total energy variation can be either positive or negative depending on the 
parameters of the model. This is different from the situation considered in [11], where 
that variation was always non- negative. 



2 Proof of Theorem 1.1 

2.1 Generator of dynamics K rn 

We recall the definition of the so-called 'C-Liouvillean' introduced for the study of open 
systems out of equilibrium in [18], and further developped in [11, 12, 13, 26, 27, 29] (see 
also references in the latter papers). Let («/#, A#) denote the modular data associated 
with (SDT#, ^#), with # given by S,TZ or £. Then 

(J, A) = (J s ®Jn® Je, A 5 ® A TC <g> A £ ) 
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are the modular data associated with (Wts <£> Tin &> Tie, ® ® ^e)- We will 
write J m and A m to mean that these operators are considered on the m-th copy of the 
infinite tensor product Tic- 

We define the C-Liouville operator 

K m = L m - J m A^y m (A)A m 1 /2j m = L + W m (X), 

m > 1, where W m (X) G 9Jt 5 ® 9% <8> OJfc is given by 

W m {\) = X^Vsn-iJs&T^Jn^WsniJs&T^Jn^l 2 )) 

+\2(V S£ ,m - (Js&T ® J m AV 2 )^ ifn (J.sAy 2 ® JmA, 1 ,/ 2 )) 
= AiWsw + A 2 W5 £iTO . (2.1) 

Of course, Wss trn is the operator acting as Wse on the subspace TLs <S> of "H, and 
trivially on its orthogonal complement. 

The operators K m have two crucial properties [18, 11, 26]. The first one is that 
they implement the same dynamics as the L m : 

e itim Ae~ itLm = j tKm Ae~ itKn , Vt > 0, VA G Tl s ® Tin ® 9tt C - 

The second crucial property is that the reference state (1.9), is left invariant under 
the evolution e ltKm , 

K m V = 0, Vm. (2.2) 
2.2 Reduced Dynamics Operator 

We follow the strategy of [11] to reduce the problem to the study of the high powers of 
an effective dynamics operator. The main difference w.r.t. [11] is that in the present 
setup, the effective dynamics operator acts now on the infinite dimensional Hilbert 
space Hs <8> Hn- 

We first split off the free dynamics of elements not interacting with S by writing 
the product of exponentials in U(m), (1.4), as 

U(m) = U m e~' 1TLm e^ TLm -' ■ ■ ■ e~ irLl U+, 

where 

Lj = L S + L n + L e + V{\) 
acts nontrivially on the subspace TCs <8> Ti-n <8> Ti.£ j and 

m j—1 

U m = eX p(- iT ^^ L ^)> 

3=1 k=l 
m j—1 

U+ = exp(-ir^^L ffe ). 

3=1 k>j 
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Let O m be an instantaneous observable (1.7). A straightforward computation shows 
that (see also [13], equation (2.19)) 

a m (O m ) = (C/+)*e iTil • • -j TL ™N(O m )e- WL ™ ■ ■■e~ iTLl U+, (2.3) 

with 

M{O m ) = Ask ^m+j^^B^^) ® r j=0 # m+j ( B j)- ( 2 - 4 ) 

As normal states are convex combinations of vector states, it sufficient to consider 
the latters. Let V& p be the GNS vector representing the initial state p, i.e., p(-) = 
{Wpl • typ). It follows from the separating property of (see (1.9)) that given any 
e > 0, there is a B' G 9Jf s.t. \\^> p - B' ^> $\\ < e. Next, we approximate B' = B' + b e (N), 
where b e (N) — > as N — > oo (for each e fixed), and where B' G 971' has the form 

B' = B' S ® B' n ®% =1 # n ® n > at +1 l £n , (2.5) 

with 5^ G 9Jt'_^. The vector -B'^o is thus approximating the initial state p . Let O m 
be an instantaneous observable and let us consider the expression 

(B'y \a m (O m )B'y ) = (* \(ByB'a m (O m )y )- 

We use expression (2.3) and the properties of the generators K n to obtain 

(V \(B')*B'a m (O m )V ) = (^ \(B')*B'(U+r^ TKl ■ ■ ■J rKm tf(O m )V ). 

Note that U+^> = *o- Let 

P N = l s ® l n ® %! ® • • • 1^ ® ^* £jv+1 ® ^ £jv+2 ® • • • > 

where -P* £fe = \^e k ){^e k \- Since {B')*B' acts non-trivially only on the factors of the 
chain Hilbert space having index < N, we have 

(y \(B')*B'a m (O m )V ) = (2.6) 
(*o\(B')*B , (U%)*e iTKl ■ ■ ■ e irKN P n z 1tKn + 1 ■ ■ ■ e iTKm Af{O m )^ ), 

where (for m > N; we have the limit m — > oo in mind) 

N-l N 

U+ = P N U+ = exp(-irJ2 L ^)- 

j=l k=j+l 

Recall 

Psn = % ® In ® |*c)(*c|- 

We have for m > N + I 

P N j rKN + 1 ■ ■ ■ e irKm Af(O m )^ = Psu^ tKn+1 ■ ■ ■ e irX "W(O m )^ 

= PsnM™- 1 -"- 1 ^^™- 1 ■ ■ ■ e irKm M(O m )^ , (2.7) 
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where we have introduced the following reduced dynamics operator (RDO), see (2.21) 
in [13] 

Psn^ TK Psn = M <g> |* c )(*c| - M acting on H s ® Hn- (2-8) 
In the last step of (2.7), we use the property 

Psnj TKs j TKs+1 ■ ■ ■ j TKt Psn = Psne [TK °Psne iTKs+1 Psn ■ ■ ■ Psk^'Psk, 

which holds for any 1 < s < t. This property follows from the independence of the £j 
for different j, see [11], Proposition 4.1. 
Combining (2.6) with (2.7) we obtain 

(y \(B l )*B'a m (O rn )y )= (2.9) 
($o\(B')*B'(U£)*e 1TKl ■ ■ ■ e irKN P S nM m - l - N - 1 e [rKm - 1 ■ ■ ■ e iTKm Af{O m )^ ). 

In order to emphasize the dependence on the coupling constants A = (Ai,A2) we 
write K(X) and M(A). The following are general properties of the RDO. 

Proposition 2.1 Let A E M 2 be arbitrary. We have 

i) M{\)eB{H sn ) 

ii) M{\)^ sn = $> sn , where $> sn :=^ s ®ty n 

hi) For any ip in the dense set V = {Asn ^sn, ^STl £ ^sn}, there exists a constant 
C(ip) < oo s.t. 

sup||M(A)>|| < C(<p). (2.10) 

nGN 

Proof: i) follows from the fact that K is a bounded perturbation of a self-adjoint 
operator and ii) is a consequence of (2.2). To prove hi), hrst note that T> is dense since 
^57^ is cyclic for Wlsn ■ Then note that the following identity holds for all Bsn £ 9Jlsn 

(B sn ^o\a n (A sn )^ ) = (B S n^sn\M(X) n A sn ^ sn ). 

Statement hi) of the lemma follows from the density of T> and unitarity of the Heisen- 
berg evolution, with C (Ag-ji^f sn) = HAs^H- 

Remark. Contrarily to the cases dealt with in [11], [12] and [13], where the underlying 
Hilbert space is finite dimensional, we cannot conclude from (2.10) that M(A) is power 
bounded. Hence we do not know a priori that <r(M(A)) C {z : \z\ < 1}. 

2.3 Translation analyticity 

To separate the eigenvalues from the continuous spectrum, we use analytic spectral 
deformation theory acting on the (radial) variable s of the reservoir 1Z. 

Recall the definition (1.10) of the translation. It is not difficult to see that 

K e := T{6)~ l KT{6) = L o + 0N + \xW S n(0) + MW S£ , (2.11) 

where N is the number operator, and that the right side of (2.11) admits an analytic 
continuation into 9 £ Kg , strongly on the dense domain D(Lq) n D(N), defining a 
family of closed operators (see [17]). 
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Theorem 2.2 (Analyticity of propagator) Assume that HI and H2 hold. Then 

1. T(9)~ 1 e lrK T(9) has an analytic continuation from 6 G R into the upper strip 
kq , and this continuation is strongly continuous as Q9 J, 0. 

2. For each 9 G k$ Ul, the analytic continuation of T(9)~ l e lrK T(9) is given by 
e lrKe , which is understood as an operator-norm convergent Dyson series (with 'free 
part' e iT ( L °+ dN )). 

3. For each 9 G k$ , \j i— > T(9)~ 1 e 11 ~ K T(9), j = 1,2, are analytic entire functions. 
Remarks. 1. The proof of this result yields the following bound for all 9 G ng U R, 

\\T(9)~ 1 e irK T(9)\\ < e rsup o<ae<e () ll^ll. 

2. If 0i,02 £ K8o W1 th 0i + 02 G Ke , then 

T(0i + 2 )~ 1 e irA 'T(0i + 2 ) = T(0 2 )- 1 T(0i)^ 1 e iTA 'T(0i)T(0 2 ). 

In particular, T(9)~ 1 e lrK T(9) is unitarily equivalent to T(iQ9)~ 1 e lrK T(iQ9), via the 
unitary T(K0)- 1 . 

Proof of Theorem 2.2. For G R, the Dyson series expansion of T(9)~ 1 e lrK T(9) 
is given by 

T(0)- 1 e irA 'T(0) (2.12) 
= Vi"/ dti • • • / ~ di n e^^We^)^^- 1 -*™)^^^!)^^- 2 -*™- 1 )^ • • • 

n=0 ^ ^ 

... e i(t 1 -t 3 )9JV Wj?(tl)e i(r-t 1 )tfJV e lr£o j 

where we define 

W" fl (t) = e itio W e- itLo 
In the derivation of (2.12), we use that for all t G R, 

T(9)~ 1 e itLo T(9) = e it ( L o+ eN ) = e itLo e it9N _ e it6N e itL ^ 

All the operators e 1 ^ fe ~ 1_ifc ^ Ar , as well as e ltn8N and e^ 7 " - *!)^ appearing in the integrand 
of (2.12) have analytic extensions from real to G Kg which are continuous at R, and 
each of those extensions is bounded, having, in fact, norm one (uniformly in G Kg UR 
and in the tj). Consequently, due to assumption H2, for fixed values of t±, . . . , t n , the 
integrand in (2.12) has an analytic extension into ng which is again continuous on R. 
Let's call this extension ht lt ...j n (0) ■ For G kq , we have 

\\9eh tl ,..., tn (0)\\ < C(\\dgW \\ + [90]- 1 ), (2.13) 

uniformly in t\,...,t n , for some constant C (which depends on n). It thus follows, 
using the Lebesgue Dominated Convergence Theorem, that the integral on the r.h.s. 
of (2.12) is analytic in G k$ . To show (2.13) we note that the derivative dgh is a 
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sum of terms where dg is either applied to one of the W$ or to an exponential. In the 
latter case, we have 

\\d e j tm \\ = ||tiVe-* 3eAr || = [^ri^WVe-'^H < {^e^supxe-*, 

x>0 

where dge ltdN is understood in the strong sense. 

The norm of the integral on the r.h.s. of (2.12) is bounded above by 

[t sup < Qe<6)() 

uniformly in 9 G Kg. It follows that the series (2.12) converges uniformly (Weierstrass 
M-test) and therefore the r.h.s. of (2.12) is an analytic function in 9 G ng . 

Using the Lebesgue Dominated Convergence Theorem and the Weierstrass M-test 
as above, one readily shows that the series in (2.12) is strongly continuous as 9$ j 0. 
Since equality (2.12) holds for real 9 this means that indeed the Dyson series is an 
analytic extension of T(9)~ 1 e lrK T(9) into 9 G Kg . Note that the series is indeed the 
Dyson series of e lrKe . 

Finally, analyticity in Ai and A2 is clear from (2.12) and (2.1). ■ 

The following result follows immediately from Theorem 2.2 and definition (2.8). 

Corollary 2.3 Recall the definition (2.8) o/M(A). If assumption H2 is satisfied, then 
T(9)~ 1 M(X)T(9) has an analytic continuation into 9 G Kg , denoted Mg(X), and this 
continuation is continuous at M. We have Mg(A)^sn = ^fsn- 

2.4 Convergence to asymptotic state 

In order to make a link with the dynamics of observables, we insert 11 = T(9)T(9)~ 1 
(with 9 G M) into equation (2.9) to obtain 



(y \(B')*B'a m (O m )V ) = (2.14) 
(^o\{B')*B'{U^)*e iTKl • • • e iTKN P sn T{9)M 1 9 n - l - N - 1 T(9)- 1 e iTK ^- 1 ■ ■ ■ e iTX "W(O m )^ ) • 



If O m is analytic (see definition (1.11)) then so is M(O m ) (see (2.4)), and therefore, by 
Theorem 2.2, e lrKm -' ■ ■ ■ e 1TKm Af(O m ) is analytic as well (and continuous on the real 
axis). The r.h.s. of (2.14) is thus an analytic function in 9 G Kg a (c.f. Theorem 2.2 and 
Corollary 2.3). Moreover, this function is continuous on M, and, by unitarity of T{9) 
for real 9, it is constant for real 9. Hence (2.14) is valid for all 9 G Kg U WL. 

One expects that the operator Mg{\) m ~ l ~ N ~ l converges to the projection onto the 
manifold of its fixed-points, as m — > 00. Under certain (physically reasonable) condi- 
tions, this projection has rank one and is given by \ 1 $>sn){' l PgW\i where tpg(X) is the 
unique invariant vector of the adjoint operator [M^A)]*, normalized as sizing W) 
1 . 

Lemma 2.4 Assume that Condition FGR holds (see before (1.12)). Then for < 
|A| < A , 

hm M ei (X) n = P h M ei (x) = |^)(^(A)|, 
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where denotes the spectral projector of the operator M corresponding to the eigen- 
value 1. The convergence is in operator norm, and occurs exponentially quickly. 



Proof. Using FGR we have 

M = P 1>m + Mq, where M Q = QMQ, Q = 1-Pi,m, 

with 

Pi,m = \*sk)W*\, (V*|*S7e) = I- 
Moreover, there exists 7 = 7e 1 (A) > s.t. 

ct(Mq) C {z £ C : \z\ < e" 7 }. (2.15) 
Therefore, the spectral radius of Mq satisfies 

spr(A/ Q ) = lim IIM^II 1 /™ < e~ 7 < 1. 
This, together with the identity 

M n = P 1>M + Mq 

yields for n large enough, 

||M n -Pi iA /|| <e~ n7 ^0 as n -> 00. 
(For some slightly smaller 7 than in (2.15) above.) 

It is now apparent from (2.14) and Lemma 2.4 how to complete the proof of Theorem 
1.1: the increasing power of Mq drives the system to an asymptotic state. Some care 
has to be exercised in the implementation of the complex deformation in the remaining 
part of the proof. Here are the details. 

The approximation of the initial state p by the vector -B'^O) explained in (2.5), 
gives us the estimate 

p(a m (O m )) = (* |(S / )*5 / (C/+)*e Ml • • • j TL ™N(O m )e- iTL ™ ■ ■ ■ e- irL ^ ) + Ri(O m ), 

where R\ is a bounded linear functional on 9JT with liniTv^oo 1 1 Pi II < e (uniformly in m; 
note also that ||O m || is independent of m). Here, e is the approximation parameter fixed 
after (2.4). Let a > and define the spectral cutoff operator \a '■= x(|dr(9 s )| < a), 
acting (non-trivially only) on Ti-jz- (Here, x(\ x \ — °~) equals one if |x| < a and zero 
otherwise.) The role of x& 1S to smoothen the deformation operators: indeed, Xo-T(9)~ 1 
is analytic entire in 9, see also (1.10). As a — > 00, Xa approaches the identity in the 
strong operator topology. Consequently, we have (see also (2.9)) 

p(a m (O m )) 

= (y \(B')*B'(U+)*Xaj TLl • ■ • e iTL ™AA(O m )e" irLm • • • e" iTLl * ) + R 2 (O m ) 
= (V \(B')*B'(U£)*X*j TKl ■ ■ ■ e iTKN M m - l - N - 1 e irK ^~' ■ ■ ■ e lrKm N{O m )^ ) 
+P 2 (O m ), (2.16) 
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where R2 is a bounded linear functional on 9Jt satisfying \\m. a ^ 00 R2 = Ri, and in 
particular, lim./v->oo lim^oo H-R2II < e (uniformly in m). 

We now introduce the spectral deformation in the main term on the r.h.s. of (2.16). 
We have for all 9 G R 

^> q \{B'YB'{U+YxJ tKi ■ ■ ■ e iTKN M m - l - N - 1 e iTK ^- 1 ■ ■ ■ e iTi ^AA(O m )^ ) 
= (^o\(B'yB'{U+)* X aT{6)e iTK ^ ■ ■ ■ ^K N , e M m-l-N-i x 

*T{6)- l e 1TK ™- 1 ■ ■■e iTKm Ar{O m )^ ), (2.17) 

where K jt g = T{9)~ l K j T{9). S ince Ora is an analytic observable, and according to 
Theorem 2.2, the right side of (2.17) has an analytic continuation into 6 6 Kg and 
this continuation is continuous at R. Moreover, on R, this continuation is a constant 
function (equal to the left side of (2.17)). It follows that the analytic continuation is 
constant on the whole region of analyticity plus the real axis, and (2.17) holds for all 
6 G Kg U R. Due to the Condition FGR and Lemma 2.4, we have 

(y \{B')*B'(U+yxae 1TKl ■ ■ ■ ^ K N M ^-i-N-i^rK m ^ . . . e i ^AA(O m )^ ) 
= ^ Q \{B')*B'{U+y X *T{e x )j TK ^ ... e iT *".«i* ) x 

x(r ei \T(9 1 )- 1 e iTK ^---e iTK -Af(O m )^o)+e- m ^R3(O m ) 

= ||B , * || 2 (^JT^ 1 )-V T ^---e ir ^(O m )*o)+e- m ^ 3 (Om), (2-18) 
where 7 > (see proof of Lemma 2.4) and where 

\\R3(O m )\\ < C(e ,N)\\T(e 1 r 1 e iTK ™->---e iTK ™N(O m )*o\\ e 9 ^e N \ 

The latter quantity is bounded uniformly in m. To arrive at the second line in (2.18) 
we made use of the fact that ^0 is invariant under the action of all of Xa, T{9) and 
e irK j: e "Yy e combine estimates (2.16), (2.17) and (2.18) to arrive at 

limsup|p(a m (O m )) - (re 1 \T(9 1 )- 1 e iTK ™~>---e iTK ™N(O m )* )\ 

< |||B^o|| 2 -l| S np\(r ei \T(9 1 )- 1 ^ K ^---^ K ^(O m )^ )\ 

m 

+ lim sup I R 2 {O m ) I . (2.19) 

m— >oo 

By taking in (2.19) first a — > 00 and then N — > 00, we see that the r.h.s. of (2.19) is 
bounded above by e (note that in this double limit, H-B^oll tends to \\^ p \\ = 1). Since 
e was chosen arbitrarily small to begin with, it follows that the left hand side of (2.19) 
is zero, and thus (1.12) follows. This completes the proof of Theorem 1.1. ■ 

3 Analysis of M(A) 

An important issue in the analysis of concrete models is the verification of the Fermi 
Golden Rule assumption FGR (see before Theorem 1.1). We have introduced the 
description of the two types of systems, 'small' and 'reservoir' in Section 1.1. For a 
more detailed analysis, we need to complement that description. 
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We denote eigenvalues and associated eigenvectors of Hs by 

Ei,--- ,E d , and <pi,- mm > <A*- 

Before analyzing the spectrum of Mg(\) in general, we mention some easier special 
cases. 

• In the unperturbed case (A = 0) we have 

M(0) = e irLs <g> e lrLlz , with 
ff(M(0)) = U< E i- E ")\ ufe", leRl, 

where the eigenvalues Q 1T ( E i~ E *) are embedded and have corresponding eigenvectors 
<Pj <8> (fk <H> V&ft. The eigenvalue 1 is at least d-fold degenerate. 

• In case the coupling Ai between the small system and the reservoir is zero, we 
have 

M(0, A 2 ) = M(A 2 ) ® e iTL « on W 5 ® 

where 

M(A 2 ) ^ P iS e iT ( Ls+i£+A2ys£ )P lS and P 5 = Is <8> |*£>(*fi|. 

The results of [11] apply to M(A 2 ), which is nothing but the RDO corresponding to the 
repeated interaction quantum system formed by S and C only. In particular, we get 
that M(0, A 2 ) is power bounded, as M(A 2 ) is and e lTljR is unitary. Moreover, assuming 
the interaction Vss "effectively" couples S and C, hypothesis (E) in [11], we know that 
the spectrum of M(A 2 ) satisfies 

<t(M(A 2 )) = {fij(X 2 )} j=i,2,-..,cP, 

with jUi(A 2 ) = 1 a simple eigenvalue with eigenvector tyg and //j(A 2 ) G {z \ \z\ < 1}. 
Hence, 

a(M(0, A 2 )) = a(M(A 2 )) U {|^(A 2 )|e", I G M} i= i,..., d 2, 

where the eigenvalues are embedded in the absolutely continuous spectrum again. 

• In case the chain is decoupled, i.e. if A 2 = 0, we get 

M(Ai,0) = e iT ( L s+Ln+XiW sn ) Qn -Hs^Hn, 

whose spectral analysis already requires the tools we will use for general A. 

We now turn to a perturbative analysis of Mg(X) (small A). Take 8 G Kg and let 
Ai = A 2 = 0. Then 

Mfl(0) = e^Ls+Ln+dN) = e irL s ^ jrLnjrON 

and 

a(Mg(0)) = },.,.<,.....,;; U {eV^, I G R} jeN .. 

The effect of the analytic translation is to push the continuous spectrum of Mg(0) onto 
circles with radii e~ T ^ s , j = 1,2, . . ., centered at the origin. Hence the discrete spec- 
trum of Mg(0), lying on the unit circle, is separated from the continuous spectrum by 
a distance 1 — e~ T ^ e . Analytic perturbation theory in the parameters Ai,A 2 guaran- 
tees that the discrete and continuous spectra stay separated for small coupling. The 
following result quantifies this. 
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Proposition 3.1 Let Co (A) := supg 6KfloUK ||We||. Take 6 G Kg and suppose that 

rCo(A)e rC °( A ) < ±(1 - e" r3e ). (3.1) 
Then the spectrum of Mq(X) splits into two disjoint parts, 

a(M e (X)) = af\\) U <rjp(\) with a^\\) n ^(A) = 0. 
These parts are localized as follows 



o- 



(0) (A) C {z : 1- i(l-e- T3e ) < \z\ < 1} (3.2) 



4 1} (A) C {z : 0< |z| <e- r30 + ^(l-e- r ^)}. (3.3) 

Moreover, the spectrum o~q (A) is purely discrete, consisting of d 2 eigenvalues ( counted 
including algebraic multiplicities) . 

Proof. According to the Dyson series expansion (2.12), we have 

M e {\) = e iT( - Lo+eN 1 + S 2 , (3.4) 

where 



1 



S 2 \\< V < rC (A)e T 



n=l 

Since e ir( - Lo+0N ^ is a normal operator, the spectrum of the perturbed operator Mg(X), 
(3.4), lies inside a set whose distance to a(e 1T ^ Lo+eN ^) does not exceed HS^H- The 
spectrum of e 1T ( L o+SN) cons [ s i s Q f isolated eigenvalues lying on the unit circle and 
of continuous spectrum lying on concentric circles centered at the origin, with radii 
e"™ 39 , n = 1,2, ... 

It follows that if (3.1) is satisfied, then the continuous spectrum of Mg(X) is located 
as in (3.3). Furthermore, the spectral radius of Mg(X) is determined by discrete eigen- 
values only, and so by Lemma A. 2 below, these eigenvalues cannot lie outside the unit 
circle, from which (3.2) follows. ■ 

Remark. We always have 1 € eri (A), with eigenvector ^sn- 

As a consequence of Proposition 3.1, a verification of FGR for concrete models, like 
the one of Section 1.4, is done via (perturbative) analysis only of the discrete eigenvalues 
of M e (X). 

4 Energy fluxes, entropy production 

We use the notation and definitions of Section 1.3 and assume thoughout that Assump- 
tion H3 is satisfied. We have 

a m+l {L n ) - a m {L n ) = a™«+i(^) - L n ), 
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with 

<{■) = e irZm • e" iTi "\ L m = L + # m (V). 

Thus 

a T m+l {L n ) - Ltz = i / a t m+1 ([L m+1 ,LT Z \)dt = i «L+i([^iKsRi Ln]) dt, 
Jo Jo 

and we arrive at the expression for the variation of energy in the reservoir 

AE n (m) = a m AjT a^^Vsn, Ln]) dt\ • 

From definition (1.7) and assumption H3 (see before (1.14)) together with the analyt- 
icity of the dynamics (Theorem 2.2), it is clear that o4j +1 ([AiVs7£, Ln\) is an analytic 
instantaneous observable. By approximating the integral J Q T dt by a Riemann sum 
(converging uniformly in m) we see that the integral in question is the limit of a sum 
of instantaneous observables all having uniformly bounded indices I, r (see (1.7)). Fur- 
thermore, a m is bounded uniformly in m and therefore we can apply Theorem 1.1 to 
conclude that 

Yim o p{AE n (m)) = p+, A (i J a^X^Vsn, L n }dt\ . (4.1) 

Next we examine the variation of energy in the chain C. During the time interval 
[tot, (to + 1)t), the energy of the element k ^ m + 1, of the chain is invariant. 
Hence, the variation of energy in the whole chain between the times m and to + 1 
coincides with that of the element £ m +i only. Thus, 

AE c (m) = a m +\L £m+1 ) - L £m+1 = a m +\L £m+1 ) - a m (L £m+1 ). 

Proceeding as above, we arrive at 

AE c (m) = a m fijf^^M^i^),^])^ , 

where i J Q T a t m+1 ([S m+ i(Vs £ ),L £m+1 ]) dt is an instantaneous observable (Assumption 
H3). 

Finally, for the variation of energy in the small system S we get 

AE s {m) = a m (i £ a^ +1 ([A 2 ^ m+1 (y 5 £) + L s ]) dt) . 

It follows now from Theorem 1.1 (by the same reasoning leading to (4.1), see also 
[11]) 

dE™ = P+ ,xU tot ) = p+,x(y - « r (^)), 



where 



j tot =V- a T {V) = -i [ T a\[L s + L n + L £ , X 2 V S£ + XiVsn}) dt. (4.2) 

Jo 
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On the other hand, by the same reasoning, dE+, for # = S,£,1Z (see (1.14)), is 
given by p +i \(j#), where 

f = i r a?([\2Vse + XiVsR„Ls])dt, (4.3) 
J o 

f = i [ T a t ([X 2 V S£ ,L £ })dt, (4.4) 
j o 

= i ra*([A!^,L^])dt. (4.5) 



o 



The following result relates the various flux observables j*. 

Proposition 4.1 We have j tot = j s + j £ + j n . Furthermore, dE+ := p + ^\(j s ) = 0. 

Proof. The relation j tot = j s + f + j n follows directly from (4.2) and (4.3)-(4.5). To 
sec that dE% = we note that since L$ is bounded we have 



1 1 N 

-(a N (L s )-L s ) = -Y,a m (L s )-a m - 1 (L s )^0, as iV ^ oo. 



m=l 



The main ingredient in the analysis of the entropy production is the following 
entropy production formula, established in [19] 

AS(m) 

= P (u(m)* [p £ £ L £ , k + (3sL s + PkLk] U(m) - p e L ^ k ~ ~ ^n) ■ 

k k 

Following the proof of Proposition 2.6 of [11], we get 

m—1 m—1 m 

AS(m) = (p n -p £ )^p(AE n (k)) + (p s -p £ )Y,P(&E s (k)) + (3 £ Y,P(&E tot (k)) 

k=l k=l k=l 

+Pe [p(^iVsn) - p(a m (\iVsn)) + p(hMVse)) - p(a m (\ 2 $ m+1 (V S£ )))] . 

Hence using Proposition 4.1, and since p(a m (Vsn)) and p(ct m ($ m+ i(Vs £ ))) are bounded 
in m, we get 

dS+ : = hm = _ p^dE* + &d£f . (4.6) 

m— »oo 771 

5 More detail on the concrete example 

We consider the model described in Section 1.4. In order to write down explicitly 
all interaction operators appearing in the C-Liouvillean, we need the explicit form of 
modular data of S, TZ and £. 

The modular data of S and the elements £ of the chain associated to the reference 
states ps and P/3 £:£ are given by 

Js(<f><8> VO = i>® <f>, A 5 = 11 C 2 ® 1 C 2. 
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The one of 1Z is given as follows (see also Theorem 3.3 of [18]) 



At?. = e~^ L 



J R = (_l)W)/2 Cof . 

Here, we have introduced the notation ft{s) = ie~^ S//2 fp{s) = fp(—a), where the bar 
indicates the complex conjugate. Furthermore, N = dr(l) is the number operator, 
C is the complex conjugation operator and T is the sign flip operator acting on / € 
®] =1 L 2 (R,&) as 

(J r f)(s 1 ,s 2 , ...,s n ) = /(-si, -s 2 , • • • , -s n ). 

An easy computation leads to the following expression for the "interaction" part 
W of the C-Liouvillean, 

W(X) = Ai^OlcaO^C/^Ca)) 

-%> ® ^ ® r(-l) (a*(/^(s)) - a(e-^ s //3 K (s))) ) 

+ A2 (as ® He 2 ® a £ 8) lc 2 + a 5 ® ^C 2 ® a £ ® He 2 

_ e ^^/2]l c2 a * ]l c2 g, a£ _ ® a s ® 1 C 2 ® a|) . 

Assumption (H4) on the form factor / ensures that assumption (H3) of Section 
2.3 is satisfied with 9q = S. We can thus apply the general results of Section 3. In 
particular, the map 9 1— > T(9)~ l M(X)T{9) has an analytic continuation in the strip Kg 
(see Corollary 2.3). We then fix some 9\ 6 Kg such that 1 — e - ' 1 " 3 ^ 1 ) > 0. For A small 
enough eq. (3.1) is therefore satisfied, so that we can verify the (FGR) hypotheses 
using perturbation theory for a finite set of eigenvalues, those four eigenvalues which 
are located in a^\x) (see (3.2)). When the coupling constants are turned off, we have 

aP(0) = a(e irLs ) = {l,e 1TEs ,e-' 1TEs } 

where the eigenvalue 1 has multiplicity 2. In order to make the computation in per- 
turbation theory as simple as possible, we will assume that these eigenvalue do not 
coincide, i.e. tE$ £ vrN. However, this assumption is certainly not necessary. 

Using a Dyson expansion for Mg(X) as in the proof of Theorem 2.2 and regular 
perturbation theory (see e.g. [16, 21]) we compute the four elements of a^\x). We 
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know that 1 always belongs to o~q (A). The other ones respectively write 



eo(A) 
e+(A) 



1 " A?r£v^||/(V^)|fc - Af^W ( I&LzM\ + 0(A 3> 



JtE s 



IT 



1 - A^-V^||/(V^)lr - A^-sinc 



s|||/(v^)ll^„ , x2 _ 2 l-sinc(r(^-£; 5 )) 



'-ds + y 2 r' 



t{E £ - E s ) 



+ 0(A 3 ), 



e_(A) 



-irEs 



1 - A?r£>/^||/(>/^)||| - A^sinc 2 ' T ^ ~ E ^ 



+i Af-PV 



*l ll/(V|s|)lll ds + A 2 T 2 1 " sinc(r(^ - E s )) 



s-E s 



t(E £ - E S ) 



+ 0(A 3 ), 



where sinc(x) = !HHgO anc j py stands for Cauchy's principal value. We thus get the 
following 

Lemma 5.1 Assume that \\f(y/Es)\\ & / and t(E £ - E s ) $ 2vrZ*, then (FGR) is 
satisfied. 

In order to compute the asymptotic state p+ t \, we compute the (unique) invariant 
vector tpg(X) of Mq{\)* (see (1.12)). Once again, standard perturbation theory shows 
that ^(A) = t/£(A) ®^n + O g (X) with 

A i7 t h + A 27 ri 
+ ^ a 2 7 ( - 2 - > + W 2) |U) ' 

A l 7th + A 27ri 

f2) (2") 

where 7 th and 7^ are defined in (1.19)-(1.20). Inserting the above expression in (1.13), 
this proves Proposition 1.4. 



A Some operator theory 

Our analysis of the spectrum of Mg(\) makes use of a translated version of (2.10), 
which replaces the powerboundedness of Me (A) in our setup. 

Lemma A.l Assume Asti and B$-ji are translation analytic in kq . Then 

sup \(BsR,*sR,\M?(\)AsR,y )\ < \\AsTi(e)\\\\BsTi(e)l 
Proof: Consider 

{B S n{0)^sn\a m {A sn {e))^ sn ) = {B sn ^> sn \T{e)M{X) m T- 1 {e)A sn ^sn) 

= {B S n^sn\M e (X) m Asn^sn) 
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We can use this property to bound the spectral radius of Mq{X) when it is deter- 
mined by discrete eigenvalues only. This means that there are finitely many eigenvalues 
ctj , j = 1, . . . , N, all of equal modulus a, such that sup{|z| : z £ a(Mg(X))} = a and 
a css (M e {\)) n {\z\ = a} = 0. 

Lemma A. 2 Assume that for some 6 £ Kg , spr(Mg(A)) is determined by discrete 
eigenvalues only. Then spr(M#(A)) = 1. 

This is an application of the following result stated in a more abstract setting. 

Proposition A. 3 Let M be a bounded operator on a Hilbert space Ji such that: 

i) there exists a dense set of vectors C C TC satisfying 

sup|(^|M»| < C(<p,iP), V^^eC, 

neN 

ii) spr (M) is determined by discrete eigenvalues only, i.e. 

a(M) n {z E C | \z\= spr (M)} C a d (M). 

Then, 

spr (M) < 1 

and the eigenvalues of modulus one, if any, are semisimple. 

Proof: Let {ay}j=i,- ,n, be the discrete eigenvalues such that \aj\ = a = spr (M) and 
let Pj and Dj be the corresponding eigenprojectors and eigennilpotents. Recall that 
[D^Pj] = and PjP k = 5 jk Pj, for all j, k E {1, • • • , N}. Setting Q = l- Ylf =1 Pj, we 
can write, by assumption ii) 

N 

M = J2 a j p j + Dj + QMQ, (A.l) 

3=1 

where 

\\{QMQ) n \\ < e Pn with 0<\na. 

Let K £ W be such that Df +1 = for all j £ {1, • • • , N} and D? £ 0, for some 
jo £ {1, • • • ,N}. If all eigennilpotents are zero, we set K = 0. Using the properties of 
the spectral decomposition (A.l), we get for any n £ N large enough 

Mn = E + E (f) «r fcjD ") + (Q M Q) n - 

j=l \ k=l ^ ' / 

Consider first the case if = 0, where all .Dj = 0. Assume that a > 1 and consider 
(p £ H such that Pj y> 7^ 0, for some jo E {!,■•• , -/V}. We define <^ = Pj ^/||Pj ^|| 
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such that M n ifo = cv^po- Now, C being dense, for any e > 0, there exists <po G C with 
\\<Po — ¥>o|| < e so that 

JV 

APft, = M> + £a£.P i (£ - (^o) + (QMQ) n (<p - p ), 

3=1 

where the norm of the last two terms is bounded by a n e ^X^j=i ll^j'll + e (P~ lna ) n 
Hence, 

(ip \M n if ) = o 1 j ({ l Po\fo) + 0(e)), with 0(e) uniform in n, 

and (tpo\(po) = 1 + 0(e). Thus the modulus of the RHS goes to infinity exponentially 
fast with n (since |cty | = a > 1), whereas the LHS should be uniformly bounded in n 
by assumption i). 

Consider now K > and let ip € TC be such that D^ip ^ 0. Assume a > 1 and set, 
as above, (fo = Pj ip/\\Pj ip\\. We have for n large enough 



where, for 1 < k < K and n large, 

Let ipo = D^po/WD^ipoW 2 , and, for any e > 0, <po, ipo in C such that \\tp~o — tpo\\ < e 
and ||V>o — ^o|| < £• Then, as n —> oo, 

x 3'o \ \K 1 "J" 



= t 1 + Of 1 /"))- 

Thus 

$o|AP£o> = (^o|M> ) + <^o - V-o|M"^ ) + (^o|M"(^o - <A))>, 
where the vector 



m^o = E f a l p i + E ft) a T kD ) 1 + (QMQT&o 

3=1 V fc=l ^ ' 

satisfies for n large enough and some constant C uniform in n, 

\\M n <p \\ < Ca n {j\ ||0o|| < Ca" (™) (1 + e), 
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and a similar estimate holds for ||(M n )*^o||- We finally get, for some constant C, 
uniform in n and e, 



Again, if a > 1, the RHS diverges as n — > oo whereas the LHS should be bounded by 



Remark: To get Lemma A. 2, from this Proposition, note that ^>siz is cyclic 
for Wis ® Win and the set of analytic observables Asn is (strongly) dense in Wis ® 
Wl-ji. Moreover, Lemma A.l shows that the dense set of analytic vectors of the form 
{As-r^sk} satisfies assumption i). Finally, as ^?s~R is invariant by Mg(X), the spectral 
radius is equal to 1. 



|(Vo|M^3o)| > |(Vo|M> )|-Cea 




II) 




ii), and the result follows. 
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